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Abstract 

We have considered the divergence structure in the brick-wall model for the sta- 
tistical mechanical entropy of a quantum field in thermal equilibrium with a black 
hole which rotates. Especially, the contribution to entropy from superradiant modes 
is carefully incorporated, leading to a result for this contribution which corrects some 
previous errors in the literature. It turns out that the previous errors were due to an 
incorrect quantization of the superradiant modes. Some of main results for the case of 
rotating BTZ black holes are that the entropy contribution from superradiant modes 
is positive rather than negative and also has a leading order divergence as that from 
nonsuperradiant modes. The total entropy, however, can still be identified with the 
Bekenstein-Hawking entropy of the rotating black hole by introducing a universal brick- 
wall cutoff. Our correct treatment of superradiant modes in the "angular-momentum 
modified canonical ensemble" also removes unnecessary introductions of regulating 
cutoff numbers as well as ill-defined expressions in the literature. 
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Since Bekenstein JIJ suggested that black holes carry an intrinsic entropy proportional 
to the surface area of the event horizon, and Hawking |2j] provided a physical basis for this 
idea by considering quantum effect, there have been various approaches to understanding the 
black hole entropy. One of them is the so-called "brick-wall model" introduced by 't Hooft 
@. He has considered a quantum gas of scalar particles propagating just outside the event 
horizon of the Schwarzschild black hole. The entropy obtained just by applying the usual 
statistical mechanical method to this system turns out to be divergent due to the infinite 
blue shift of waves at the horizon, 't Hooft, however, has shown that the leading order term 
on the entropy has the same form as the Bekenstein-Hawking formula for the black hole 
entropy by introducing a brick-wall cutoff which is a property of the horizon only and is the 
order of the Planck length. The appearance of this divergence [|J || and relationships of this 
"statistical-mechanical" entropy of quantum fields near a black hole with its entanglement 
entropy || and quantum excitations of the black hole J7J have been studied, leading a great 
deal of interest recently [§]. 

The brick-wall model originally applied to the four dimensional Schwarzschild black hole 
B has been extended to various situations. The application to the case of rotating black 



holes has also been done for scalar fields in BTZ black holes in three-dimensions 0, 10 and 



in Kerr- Newman and other rotating black holes in four-dimensions pi] , |12|| . In a background 
spacetime of rotating black holes, it is well known that scalar fields have a special class of 
mode solutions, giving superradiance. It is claimed in Ref. that the statistical-mechanical 
entropy of a scalar matter is not proportional to the "area"(z.e., the circumference in the 
three-dimensional case) of the horizon of a rotating BTZ black hole and that the divergent 
parts are not necessarily due to the existence of the horizon. Contrary to it, in Ref. |J, the 
leading divergent term on the entropy is proportional to the "area" of the horizon, and it is 
possible to introduce a universal brick-wall cutoff which makes the entropy equivalent to the 
black hole entropy. Moreover, it is claimed in Ref. || that the contribution from superradiant 
modes to entropy is negative and its divergence is in a subleading order compared to that 
from nonsuperradiant modes. On the other hand, for the case of Kerr black holes in Ref. [ Yl\ . 
the divergence is in the leading order but the entropy contribution is still negative. 

One may expect that the leading contribution to the entropy comes from the region 
very near the horizon as in the case of nonrotating black holes. Since the vicinity of a 
rotating horizon can also be approximated by the Rindler metric, it is seemingly that the 
essential feature of the leading contribution will be same as that in nonrotating cases. Our 
study in detail shows this naive expectation is indeed true. That is, we point out that 
previous erroneous results appeared in the literature were mainly related to an incorrect 
quantization of superradiant modes. For the case of rotating BTZ black holes, we have 
explicitly shown that superradiant modes also give leading order divergence to the entropy 
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as nonsuperradiant ones. Moreover, its entropy contribution is positive rather than negative 
found in Refs. 0, However, the total entropy of quantum field can still be identified 
with the Bekenstein-Hawking entropy by introducing a universal brick-wall cutoff. It also 
has been shown that the correct quantization of superradiant modes in the calculation of the 
"angular-momentum modified canonical ensemble" removes various unnecessary regulating 
cutoff numbers as well as ill-defined expressions in the literature. 

Let us consider a quantum gas of scalar particles confined in a box near the horizon of 
a stationary rotating black hole. The free scalar field satisfies the Klein-Gordon equation 
given by (□ + fi 2 )(p = with periodic boundary conditions 

<f>(r + + h) = <f>(L). (1) 

Here, r+, r++h, and L are the radial coordinates of the horizon and the inner and outer walls 
of a "spherical" box, respectively. Suppose that this boson gas is in a thermal equilibrium 
state at temperature fi^ 1 . Due to the existence of an ergoregion just outside the event 
horizon, any thermal system sitting in this region must rotate with respect to an observer at 
infinity. Accordingly, in order to obtain the appropriate grand canonical ensemble for this 
rotating thermal system, one should introduce an angular momentum reservoir as well in 
addition to a heat bath/particle reservoir characterized by temperature T = (3~ l and angular 
speed Q with respect to an observer at infinity f!3]| . All thermodynamic quantities can be 



derived by the partition function Z(/3, fl) given by 

Z(J3,SI) = Tre"^-^ 1 , (2) 

where :H: and :J: are the normal ordered Hamiltonian and angular momentum operators 
of the quantized field, respectively |14|, [13], ||. Here, we assume that particle number of the 



system is indefinite. 

As usual, by using the single-particle spectrum, one can obtain the free energy F((3, fl) 
of the system in the following form 



f3F = -\nZ = - J2^J2le~ P(£x ~ njx) ] k (3) 

A k 

/ - Ea M 1 + e - ^ - "^] for fermions, 

1 Ea H 1 ~ e-ft ex - aj *>] for bosons with e x - > 0. ^ 



where A denotes the single-particle states for the free gas in the system. E\ = :H: | 1a) 
and j\ = (1a | 'J' | 1a) are normal ordered energy and angular momentum associated with 
single-particle states A, respectively. The occupation number k = 0, 1, 2, • • • for bosonic fields 
and k = 0,1 for fermionic fields. Note that, if there exists a bosonic single-particle state with 
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its energy e\ and angular momentum j\ such that E\ — Qj\ < 0, the expression in Eq. (||) 
becomes divergent and so is ill-defined. In order to compute the free energy in Eq. ([3]), one 
must know all single-particle states and their corresponding values of E\ and j\ for a given 
system. 

As pointed out in Refs. [15, HJ, the quantization of matter fields on a stationary rotating 



axisymmetric black hole background is somewhat unusual due to superradiant modes which 
occur in the presence of an ergoregion. The mode solutions will be of the form, <f)(x) ~ 
e ~ tuJt + mi P ( because this background spacetime possesses two Killing vector fields denoted by 
d t and d^. Since the partition function in Eq. @) is defined with respect to an observer at 
infinity, the vacuum state to be defined by the standard quantization procedure should be 
natural to that observer at infinity in the far future. Thus, we expand the neutral scalar 
field in terms of a complete set of mode solutions as follows: 

fa) = £ / MVl<k + + E / n + €>™ ) (5) 

+ E / Mb^-mU^-m + b^Z-mUZ-J- (6) 

Here w out (x) describes unit outgoing flux to the future null infinity T + and zero ingoing flux 
to the horizon 7i + while u m (x) describes unit ingoing flux to 7i + and zero outgoing flux to 
T + . These mode solutions are orthonormal 

« u m , u°f m> ) = (u^ m , <V> = «_ m , u^,_ m ,) = 5(u - u')5 mm , (7) 

with respect to the Klein-Gordon inner product 



Z Jt=const. 



Note that M wm (x) ~ e - lu)t+im v an d we suppressed other quantum numbers. Modes with 
Cj = uj — Qh^ti < exhibit the so-called superradiance. Here fin is the angular speed of 
the horizon with respect to an observer at infinity. An observer at infinity would measure 
positive frequency for all modes and w™ m with Cj > 0, but measure negative frequency for 



wm 

u -u-m w ^h Cj < 0. A ZAMO |17] near the horizon, however, would see positive frequency 



waves for u m a; _ m with uj < as well as for u™ m with uj > 0. Hence, in the terminology 



of Ref. |15| , we adopt the "distant-observer viewpoint" for u (x) and the "near- horizon 
viewpoint" for u m (x). And the conventions are chosen so that they agree with viewpoints. 

Now the mode solutions for particles confined in the near-horizon box would be con- 
structed by linearly superpose u m and u out above as follows: 

f <,m + <W<L fo r £ > 0, 
<m + «™«™- m for Cj < 0, 
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with appropriate normalization factor. a um is chosen so that the modes satisfy the periodic 
boundary condition in Eq. ([[]). Thus, only some discrete (real) values of u will be allowed |[19|| . 
(fiu)m(x) are understood to be cut off everywhere outside the box. Note that (j) um (x) ~ 

e -icjt+im V for aU 

The inner product of these modes becomes 

f L 

(0u,m , 4>L>>m>) = Sujuj'Smm' (u - fi ™) |</w| 2 N _1 d£, (10) 

Jr + +h 

where we have used rf£ M = n^d'E and the unit normal to a t — const, surface = iV _1 (<9t + 
^o^) M - Here fi (r) is the angular speed of ZAMO's UTTf . Since fi (r) < Q H = f2o(r — r +), 



the norm of a mode solution with to > is positive if Q = uj — f2#m > 0. When uj < 0, the 
norm could be either positive or negative depending on the radial behavior of the solution. If 
the norm of 4> um (x) is negative, we can easily see that <fi_u-m{ x ) ~ e wt-*mtp ^gg the positive 
norm. Let us define a set SR consisting of mode solutions wm with to > whose norms are 
negative. Then, the quantized field inside the box can be expanded in terms of orthonormal 
mode solutions as 

i x ) + aLCm^)] + E [a-w-m^-u-mW + aL a; „ m 0!_ w _ m (a;)] , (11) 

A^SR AeSR 

where the single-particle states are labeled by A = (uj, m). The Hamiltonian operator in the 
reference frame of a distant observer at infinity becomes then 

H ^ Uj{o, u)m o) u]m -\- CL^ jrn CL LU m) ^ ] ( ^) (fl— ui— m®— uj— m ^— u>— m®— m) 

A^SR AeSR 

= £ 0,(i\U + i) + E (-^)(^-m + h, (12) 
A^SR Z AgSR Z 

where A^ m = a^ m a wm and iV_ w _ m = a[ u _ m a_ [d _ m are number operators. Now, by following 



the standard procedure for defining a vacuum state and single-particle states [16, [15| , we can 



easily see that (e\,j\) = (u,m) for single-particle states A = (uj,m) SR while (e\,j\) 



(—uj, —m), instead of (o;,m), for single-particle states A = (cu,m) 6 SR. In Refs. || |TI| |I~2 
however, (e\, j\) = (u,m) for A G SR have been used, and our study shows that this error 
comes from the incorrect quantization of superradiant modes. This important difference 
is a peculiar feature of the quantization of matter fields in the presence of an ergoregion 
and turns out to make our "angular-momentum modified canonical ensemble" in Eq. (0) 
being well-defined as shall be shown below in detail. It also makes somewhat unphysical 
treatment of superradiant modes and introduction of various cutoff numbers unnecessary in 



the calculation of statistical- mechanical entropy in the literature. For example, in Ref. [10|, a 
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cutoff in the occupation number k was introduced to avoid the divergent sum for superradiant 
modes in the log in Eq. (H). 

In general, uj is discrete due to the finite size of the box, but the gap between adjacent 
values goes small as the size of the thermal box becomes large. In this continuous limit, one 
may introduce the density function defined by g(uj,m) = dn(uj , m) / duj where n(uJ,m) is the 
number of mode solutions whose frequency or energy is below uj for a given value of angular 
momentum m. Thus, g(uj, m)duj represents the number of single-particle states whose energy 
lies between uj and uj + duj, and whose angular momentum is m. Using this density function, 
the free energy in Eq. (§) can be re-expressed as 

(3F = -J2[dco g(uj,m)\nJ2[e' H£x ' ajx) ] k . (13) 

The angular speed Q in Eq. ( |I"3l) is a thermodynamic parameter defined, in principle, by 
its appearance in the thermodynamic first law for the reservoir, namely TdS = dE—QdJ+- ■ ■. 
Since a particle cannot move faster than the speed of light, its angular velocity with respect 
to an observer at rest at infinity should be restricted. The possible maximum and minimum 
angular speeds are 



Ot(r) = n (r) ± yj(d t ■ djd~ v ■ d v ) 2 - d t ■ d t /d v ■ (14) 

respectively. We see that, as r — > r+, the range of angular velocities a particle can take on 
narrows down(z.e., £l±(r) — > £Ih), and so the angular speed of particles near the horizon will 
be Qh- For a rotating body in flat spacetime, one knows that all subsystems must rotate 



uniformly when the body is in a thermal equilibrium state [|13| . In fact, it is a part of the 
thermodynamic zeroth law. In curved spacetimes, the uniform rotation of all subsystems in 
thermal equilibrium may not be true to hold any more. However, since we will be finally 
interested in the quantum gas only in the vicinity of the horizon, we shall assume Q = Qh 
below. 

Now, one can see that the sum in the log in Eq. ( |i~3"D is defined well for states belonging 
to SR since E\ — Qj\ = — {uj — flam) > by the definition of the set SR. For some states with 
u — VLhTti < not belonging to SR, however, the sum becomes divergent. As will be shown 
below explicitly, however, the main contribution to the entropy of the system comes from 
the infinite piling up of waves at the horizon. For such localized solutions near the horizon, 
the signature of norms in Eq. ([H]) will be determined by that of uj — VLnm since fio( r ) — &h 
for the range of integration giving dominant contributions. Therefore, we assume that all 
states with u — f2#m < belong to the set SR. Then, the free energy in Eq. fll3"D can be 
written as F = F NS + F S r. Here 

PF NS = W 'dujg(uj,m)Hl-e-^ UJ - n « m \ (15) 
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PF SK = ^[dug(u,m)Hl-e P{U - nBm) ]- (16) 
xeSK 

In general cases, it is highly nontrivial to compute g(u,m) exactly except for some cases 
in two-dimensional black holes |p0[] . For suitable conditions, however, one can approximately 
obtain g(u, m) by using the WKB method as in the brick-wall model ||. For simplicity, let 
us consider a scalar field in a rotating BTZ black hole in 3-dimensions |[21|| . For the case 
of Kerr black holes in 4-dimensions, although the essential result is the same as that in the 
case of BTZ black holes, it requires some modified formulation basically due to the fact that 
the geometrical property near the horizon changes along the polar angle . The metric of 
a rotating BTZ black hole is given by 

ds 2 = -N 2 dt 2 + N~ 2 dr 2 + r 2 {d(f - Q dt) 2 , (17) 

where 

N 2 = r 2 /l 2 -M + J 2 /4r 2 = (r 2 - r\){r 2 - r 2 _)/r 2 l 2 , (18) 

and the angular speed of ZAMO's is Qq = J/2r 2 . Here r± denote the outer and inner 
horizons, respectively. Note that Qh — ^o( r — r +) — r_/r + Z. Then, mode solutions are 
(fiumix) = (pujm( r )^~ luJt+im{p ■ Here the radial part (p^mir) satisfies 



rN 2 — 
dr 



rN 2 ^ m (r)] + r 2 N A k 2 (r; u, m)0 wm (r) = 0, (19) 



where 

k 2 (r; u, m) = N'%oo - VL + m)(u - fi_m) - /j 2 N 2 }. (20) 

Here 0±(r) = Qo(r) ± N/r for a rotating BTZ black hole in Eq (|14D. 

In the WKB approximation, the discrete value of energy u in Eq. ( |19"D is related to 
n(uj,m) as follows 

7m(u;,m)= / dr "fc" (r;u,m), (21) 

Jr + +h 

where "fc" (r;uj,m) is set to be zero if k 2 (r;uj,m) becomes negative for given (u, m) [[|. 
Since "fc" (r; uj, m) ~ N~ 2 and N(r) — > as one approaches to the horizon, we can easily see 
that the dominant contribution in Eq. ([H]) comes from the integration in the vicinity of the 
horizon as the inner brick- wall approaches to it( i.e., h — > 0). 
Now, Eq. fll5|) becomes 



PFm = Y, fairl- T dr"k"(r;u,m) 
C^J oco Lvr Jr++h 



1 r L 

71 Jr++h 



(3 f L , [j k(r;cu,m 



ln[l - e 



-(3(UJ-n H m) 



i dr^ Iduj ■ „ , . , 
7T J r+ +h j^J eP(w-n B m) - i 

+1 ( L dr V Mr- lo, m) lnfl - e^ u - Q ^} (22) 

71 Jr + +h V W min (m) V ^ 
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by using the integration by parts in oj. For convenience, one can divide Fns into two parts 

NS - -TNS + -^NS > l^J 

where 



-,(m>0) 
NS 



l r L 



dr N 



-2 



oo /-oo 



7T Jr++h 

from states with positive angular momenta and 



dm/ du 



\J {oj — n + m)[u — c_m) 

e P{<jJ-Q, H m) _ i 



(24) 



->(m<0) 
NS 



1 



dm / oL; 



+ 



\J {uj — n + m)(u — n_m) 



7T \ Jr++h 

C ^ ^ /I ^ ^ ""J — ^ ) - 1 
— ^ drN~ 2 [° dm Jfl + fl_m 2 \n[l - e ^ Hm ) 

7Tp Jr++h J~oo V 



(25) 



from states with negative angular momenta. r erg = yMl is the radius of the outer boundary 
of the ergoregion where f2_(r = r erg ) = 0. Here we considered a massless scalar field for 
simplicity Similarly, from states belonging to SR, Eq. flTBJ ) becomes 

i 

dm 



SR 



-2 



drN 

7T Jr + +h 

+— / drN~ 2 

Tip Jr + +h 



n-m [u — n + m)(u — n_m) 

^ e -/3(0J-Q H m) _ J 

/•oo i 

jf dm^jQ+Q^m 2 ln(l - e -^ m ). 



(26) 



Note that g = —dn/duj for A G SR, and that the boundary term in Fsr exactly cancels that 



m ^ NS . 



From Eqs. (plj-Pq), we can obtain leading order dependence on the brick- wall cutoff h for 
the free energy as follows; 



7 (m>0) 
NS 



7 (m<0) 
NS 



SR 



C(3) rjl* 
f3 3 [rl-rlf 

C(3) rjP 
(3 3 {r\-rl) 2 

C(3) r\l* 

p3 ( r 2_ r 2)2 



r+ - r 2 [f— r _ r+ 
,-— In — 
h 7i n 



2V2 



$[Vh) 



n — r_ 



1„(^) + ^) 



r 2 _ r 2 ^ 



2v/2 



^2 ^2 



2vrr_ V /i y 



(27) 
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The entropy of this boson gas which is assumed to be in thermal equilibrium with the rotating 
black hole can be obtained from the free energy by using the thermodynamic relation, S = 

(3 2 dF/d(3\^ H = SPF\^ H ; 



S 



NS 



S: 



303) 
4tt 2 / 

303) 
4tt 2 / 



^,2 



- r 



2V2 



2V2 



7T h 



1 /' + 

— In — 

7r n 



ln( 



h 



ln( 



+ ^(v^) 



where the temperature of a rotating BTZ black hole |2l 

ft = (rl - rt)/27rr + l 2 . 
Now the total entropy of the system becomes 

303) 



is 



S 



r 2 ff 



(28) 



(29) 



(30) 



4V2?r 2 / V h 

In Ref. 0, it is claimed that the contribution from superradiant modes is a subleading 
order compared with that from nonsuperradiant modes. In our results above, however, we 
find that superradiant modes also give a leading order contribution which is in fact exactly 
same as that from nonsuperradiant modes in the leading order of Jr+jh. It should be 
pointed out that the entropy associated with superradiant modes is positive in our result 
whereas it is negative in Refs. 0, [T2]]. In addition, since the log terms in Eq. ( fZ8D are exactly 
cancelled, our result for the entropy of quantum field smoothly reproduces the correct result 
in the non-rotating limit(ie. ; J — > or r_ — > 0) whereas the entropy obtained in Refs. ||. |i~2 
becomes divergent in that limit. 

If we rewrite the entropy in terms of the brick-wall cutoff in proper length defined as 
^ = Ir+ +k y/9n-dr, Eq. (|30|) becomes 



S 



3Q3)/8tt 
h 



C + d{h), 



(31) 



where C = 2nr + is the circumference of the horizon. Thus, by recovering the dimension and 
introducing an appropriate brick-wall cutoff 

303); 



h 



167T 3 



l P ~ 7.3 x lQ- a l P 



(32) 



which is a universal constant, one can make the entropy of quantum field finite and being 



equivalent to the Bekenstein-Hawking entropy of a rotating BTZ black hole |21 

S = Anr + /l P = S BH 



(33) 
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in leading order. Here lp is the Planck length. For a fermionic field, although modes with 
uj < do not reveal superradiance, it turns out that only the overall numerical factor in 
Eq. (^y) is different. As mentioned before, the extension of our study to the case of Kerr 
black holes in four- dimensions is straightforward, but requires some modifications mainly due 
to the polar angle dependence of the near horizon geometry. A calculation in the phase space 
shows that the essential feature of the leading order divergence in the entropy of quantum 
fields is same as that of the present case p2| . 

Other thermodynamic quantities of quantum field such as the angular momentum and 
internal energy can also be obtained as follows; 



= 3«W 2^ + 

p=p H ,n=n H h I 



_dF 

■^matter 

Here the derivative with respect to Q has been taken for Eqs. If we put the cutoff 

value in Eq. (0), we have 

•^matter = "7 = JbR- (35) 

The internal energy of the system with respect to an observer at infinity is 

-\- QjJ J matter 



E - w m 



3^3)7167^ 4 Tj + lrg. - 4 2r 2 

- 1 3^^ + ^ ln/i )-3 M BH"3^ (36) 

where the black hole mass is Mbh — M — (r 2 + + r 2 J)/l 2 . One can easily see that J ma tter — ► 
and E — > |Mbh in the limit of non-rotating black holes (e.g., Jbh = </ — > 0). Therefore, we 
find that the entropy and angular momentum of quantum field can be identified with those 
of the rotating black hole by introducing a universal brick-wall cutoff although the internal 
energy is not proportional to the black hole mass. 

What kind of relationships could be held among parameters characterizing a rotating 
black hole and thermodynamic quantities of the system of quantum fields in equilibrium 
with the black hole? To see this, let us consider a system whose free energy depends on 



the temperature such that F((3, Vt, M, J) = (3 f(Q, M, J) f23|. Suppose the entropy of the 



system is identified with that of the black hole after an appropriate regularization, 

S = (3 2 (%) = S BU = 4vrr + . (37) 

The internal energy of the system with respect to a "corotating" observer is 

2 S Arl-r 2 _ 



E' 



p=i3 H ,n=n H 3f3 H 3 / 2 



(38) 
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Now suppose the angular momentum of the system is proportional to that of the rotating 
black hole, 



^matter i . ,, , i 

\ Oil J 



= aJ = a — ^ — . (39) 



The internal energy of the system with respect to an observer at infinity is then 

E = E' + fi^Jmatter = -— " J 2 — -, (40) 

which becomes proportional to the mass of the black hole, M = {r 2 + + r 2 _) /I 2 , only if a = 4/3. 
Therefore, we expect the relationships are probably 

■/matter = E= \ M ' ^ 

If we apply the same argument to the case of Kerr black holes in 4-dimensions, we obtain 

■/matter = J J, E =\m (42) 

of which the second relationship has been explicitly shown for the Schwarzschild black hole 
in the brick- wall model by 't Hooft ||. 



We have not obtained the relationships in Eq. ( [41]) at the present letter. The reason for 
these discrepancies is not understood at the present. It will be very interesting to see how 
the Pauli-Villars regularization method, which does not require the presence of a brick-wall 
as shown in Ref. for the case of a charged non-rotating black hole, works for the case of 
rotating black holes. 
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Especially, GK would like to thank J. Samuel, T. Jacobson, and H.C Kim for many helpful 
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